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Abstract 

Dynamics in delayed differential equations (DDEs) 
is a well studied problem mainly because DDEs 
arise in models in many areas of science including 
biology, physiology, population dynamics and 
engineering. The change of the nature in the 
solutions in the parameter space for a network 
of Phase-Locked Loop oscillators was studied in 
Symmetric bifurcation analysis of synchronous 
states of time-delayed coupled Phase-Locked Loop 
oscillators. Communications in Nonlinear Science 
and Numerical Simulation, Elsevier BV, 2014, 
(on-line version), where the existence of Hopf 
bifurcations for both cases, symmetry-preserving 
and symmetry-breaking synchronization was well 
stablished. In this work we continue the anal¬ 
ysis exploring the stability of periodic solutions 
emerging near Hopf bifurcations in the Eixed-point 
subspace, based on the reduction of the infinite¬ 
dimensional space onto a two-dimensional center 
manifold. Numerical simulations are presented in 
order to confirm our analitycal results. Although 
we explore network dynamics of second-order oscil¬ 
lators, results are extendable to higher order nodes. 

Keywords: Bifurcation, stability, time-delay dif¬ 
ferential equations, symmetry, oscillators network. 


1 Introduction 

We consider the Eull Phase model introduced in [2] 
to analyse stability of periodic orbits near Hopf bi¬ 
furcations emerging in the parameter space (/i, r) 
at the Eixed-point subspace, which are non degen¬ 
erative for the case K > 1. It has been shown 
that these bifurcations can cross the imaginary axis 
in both directions, from the left to the right and 
from the right to the left. The main approach 
used for the analysis is the decomposition of the 
infinite-dimentional space into a 2-dimensional cen¬ 
ter space corresponding to the imaginary critical 
simple eigenvalue A = iiw, w > 0, and an infinite- 
dimentional space “orthogonal” to the first one (the 
orthogonality condition will be defined below). We 
will follow closely the theory and procedures pre¬ 
sented in [iiiiiniiiiii]. 


2 The Full-phase model 

In [2] the Eull-phase model was used to find Hopf 
bifurcations in the parameter space (/r, r), the gen¬ 
eral model for a Wnode, fully-connected, second- 
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order oscillator network is: 


where (j^jr ■= (f'jit — t). Truncate the Taylor series 
up to the third-order term: 


Mt) + fj'Mt) - fi- (j)j) = 0 , 


N 


N -I 


3 = 1 


z = 1,..., A/", where: 


- r) - (j)i{t)) + - t) + 


The equilibria (f>^ in equation Q are: 
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+ </'2) sin20=^ 

“ \ -(t>if + {4>jr + cos 2<l)^ 
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here, for the sake of notation we changed 6(j)i —)■ ipi- 
Note that (j)i e C([—r,0),K), i = 

Choosing = (l)i and = ^i, the vector field 
form, is: 




i) 
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•(d (d , 

2^2 = -^^^2 + - 


N 


N-1 


3 9^-i- 
3 = 1 


(—1 -I- cos 2(j)'^)xi'^ 


n S Z, itr > 1. For our analysis, we consider three 
main assumptions: 

(a) The critical eigenvalue A of the linearization 
of Q at equilibria crosses the imaginary axis 
with non vanishing velocity, i.e. Re(A'((^^)) 

0 . 


(1 -I- cos2(j)^)x^^^ — ^ -I- sin2(/) 

(^x'f^ - -I- (^XiJ -I- cos2(/)=*= 


( 6 ) 


We define X := C([—r, 0], K^^), the Banach space 
of continuous functions from [—r, 0] into 
equipped with the usual norm 


(b) The purely imaginary eigenvalue A = iw is sim¬ 
ple. 


(c) The linearization of Q at equilibria, has no 
eigenvalues of the form ikuj, fcsZ — {!,—1}. 

The linerization of 0 at equilibria is: 


S(j)i + ^6(j). 
' N -1 




j^i r=l 
3 = 1 


5(j) 
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dcj)‘ 






(t)[=<t)± 
4>'. =<f>± 
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||x|| = sup |a;(6»)|, a; S C([-t, 0], K^^), 

-T<e<o 

and X = (x^^\... ,x^^'^) G X, where = 

{xi\x^^^). 

Now, in order to build the decomposition of the 
infinite-dimensional space, we need to define de ad¬ 
joint operator associated to the linear part of the 
linearization and a inner product, via a bilinear 
form. 

Following 013], we can represent the dynamics 
in 0 by the abstract differential equation: 

^xt{(l)) = A{r])xt{(l)) + (7) 
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which satisfies {T{t)<j)){0) = {xt{(j))){0) = x{t + 9), 
where T{t) is a semigroup of family of operators, 
9 G [—r, 0], and ry is a vector of parameters. The 












linear operator A{r]) G Mat(2iV) is defined in equa¬ 
tion (2.16) in [ 2 ], and 




|f(0) ,-r<0<O 

F{x{0),x{-T),r]) ,e = 0 

( 8 ) 


F = where /(*) = 

= 0 , and = 



Associated to the linear part of Q, we define the 
formal adjoint equation: 


= K{'n)y{t.v) + K{'n)y{t + (9) 


3 The Fixed Point space S^r 


Due to the SAr-symmetry of Q the space where 
solutions (pi lie can be decomposed into the Fixed- 
point subspace where symmetry-preserving solu¬ 
tions emerge and a subspace with symmetry¬ 
breaking solutions, this was shown in [2]. We an¬ 
alyze stability of the periodic solutions near Hopf 
bifurcations in the Fixed point space, these bifur¬ 
cations satisfy assumptions 0-0 for K > \. In 
this subspace, equation 0 has the form: 


Xi = X2 

X 2 = —y^X 2 + A '^(—1 + cos 2p^)xi 

+ K^{1 + cos 20=*=) a:I t- — -( 0 : 1 ,- -I- xiY sin 20=*= 

1 ^ 

- g [{xiT - Xi)*^ -I- {Xir -I- Xi)*** COS 20=*=] , 

( 12 ) 


then matrices Aq^t]) and Ar{ri) in ( 10 ) become: 


matrices ^ 0 ( 77 ), At{t]) are defined in [2]. The 
strongly continuous semigroup {T* {t)'ip){0) = 

{yt{tp)){0) = y{t + 0), defines the infinitesimal gen¬ 
erator: 

{A*{'q)iP) = 

%{0) ,O<0<T 

Ao{ri)'^-ipiO) + Ar{r]Y' iP{t) ,0 = 0 

( 10 ) 

such that ■^T*{t)'ip = A*T*{t)'ip, ip G X* := 
C([0,r],The natural inner product, follow¬ 
ing 0 . has the form: 

{x,y) = x'^{0)y{0) + J x{s + T)Ar{r])y{s)ds. 

For ip G X and tp G X* we have [3]: 

1. A is an eigenvalue of A(y) if and only if A is 
and eigenvalue of A*{r]). 


( i^Ai(-l-hcos(20=t)) -It)’ 

" ( Kfi{l + cosi2p^)) 0 ) ’ 

and F in 0 takes the form F = (/i / 2 )*^, with 
/i = 0, and 02 : 

f2{xt,v) = -^(xir +xi)2sin20=^ 

X 

- g [(Xlr - Xi)*^ -F {Xir + Xi)*^ COS 20=^] . 

(15) 

We need the complex eigenfunctions As('d) = 
icijs('0), A*n{0) = iu!n{0), associated to the crit¬ 
ical eigenvalues A = iw, and A = — iw with 
s('d) = si(i?) -I- 152 ( 1 ?) and n{0) = ni{0) + 1712 ( 0 )- 
These eigenfunctions can be computed solving the 
boundary value problem pjAsi ,2 = Tws 2 ,i(i?), and 
^ni _2 = ±ws 2 ,i(i?), which, after substituting the 
operator ^( 77 ), becomes: 


2. If ipi,...,ipd is a basis for the eigenspace 
of A(i 7 ) and -ipi,... ,ipd is a basis for the 
eigenspace of A*{r]), construct the matrices 
$ = {ipi,...pd) and = {ipi,...,ipd)- De¬ 
fine the bilinear form: 


^o(i7)si(0)-H A^(r7)si(-T) = -a;s2(0) . . 

Ao{ri)s2i0) + Ar{ri)s2i-T) = a;si( 0 ) ^ 


^0 ivhiiO) + {T])ni{-T) = 01712 ( 0 ) 

A]]’(77)112 (0)-f (77)112 (-r) = -01771 (0) ’ 

(17) 


(vl/,ci>)=/ 


( 11 ) 
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and from equation (21) 


with general solutions: 

= cos(wi?)ci — sin(wi?)c2 
S2(^^) = sin(wd)ci + cos(wi?)c2 

ni{6) = cos{ujO)di — sml(jj9)d2 

71,2(0) = sin(a;0)di + cos(w0)(i2 

The coefficients ci = [cn €12]^, C2 = 
[c2i C22]^, di = [dll di2]^, d2 = [d2i ^22]^ can be 
obtained by considering the boundary conditions, 

/ Aoir])+cos{u}T)AriT]) *^ 1^=0 

ujI + sm{uJT)Ar{f]) ) \ C2 ) 

( ^0 (0)+ cos(a;T)A^(7?) ^ ^ n 

\ -wl - sm{uJT)A'^{r]) ) \ d2 ) 

(19) 


the “orthonormality” condition (s, n) = 7, and set¬ 
ting Cn = 1 and C 21 = 0, see [3 [7] for more details. 

It is also possible to decompose the solution 
a;t(d) to equation ([ 7 ]) into xt{9) = yi{t)si{d) + 
2 / 2 (i)s 2 ('!^) + w(t)(d), where yi and j /2 He in the cen¬ 
ter subspace, such that yi, 2 {t) = («i, 2 ( 0 ),a:t( 0 )), 
and w in the infinite-dimensional component sub¬ 
space, thus, we have 


2/1 = W2/2 -f nj{0)F 

2/2 = -ujyi+n2{0)F 


( 20 ) 


w = A{r]) +F(xt,r])- nf (O)Fsi - (0)J^S2, 

( 21 ) 

where 
F = 

Jo, d e [-r,0) 

1 ^( 2 /i(i)si( 0 )-b 2 / 2 (t)s 2 ( 0 )-b w(t)(0)), d = 0 

( 22 ) 


3.1 The center manifold 


Following 013 HD], we know that w can be approx¬ 
imated by the second-order expansion: 

w(2/i,2/2)(^?) = ^ihi{i})yf + 2h2i'&)yiy2 + h3(d)yl), 

(23) 


thus, by differentiating and substituting equa¬ 
tion (21) keeping up to second order terms, we ob¬ 
tain: 


dw 

dt 

A{r])w + F{w + yisi + 2 / 2 S 2 ) - (di 2 Si + d 22 S 2 )/ 2 - 

(25) 


From the definition of A{r]), equivalent to (10), we 
see that 


A(ji)w = 

J \{hiyl + 2h2yiy2 + h3yl), dG[-T,0) 

\ 2lo(?7)w(0)-I-^r(?7)w(-T), d = 0 

(26) 


then, from equation ( |23[ ), ( |24[ ), ( |25[ ), and ( |26[ ), we 
can obtain the unknown coefficients hi, / 12 , and / 13 , 
solving: 


hi — 2{—ujh2 + fi^{di2Si{'d) + d22S2i'd))), 

h2 = Uj{hi - hs) + /2^(dl2Sl(l?) -b d22S2{^)), 
hs = 2(w/i2-b/2^(di2Sl(d) -b d22S2(0))), 

(27) 


Ao{v)hi{,0) + Ar{r])hi{-T) = 

2(— 01 / 12 ( 0 ) + /2^(dl2Sl(0) -b 1 / 2252 ( 0 ))), 

Ao{'n)h2{0) + Aoh2{-T) = 

O;(/li(0) — /l3(0)) -b /2^(dl2Sl(0) -b 1 / 2252 ( 0 ))), 

Aoiv)h3{0) + Ar{ri)h3(-T) = 

2 ( 01 / 12 ( 0 ) -b /2 ^(i/i25i( 0) -b l/2252(0))), 

(28) 


where = 


Id^f 


1(PJ 


^11 d^f 

0 dyidy2 


and /02 = 


2 92/1 _ 

Equation (27) is written as the inhomogeneus dif¬ 
ferential equation: 


dh 

M 


Ch -b p cos(oii?) -b q sin(oid) 


(29) 


w = -ujh 2 yl + oi(/ii - /i3)2/i2/2 + w/ 122/2 + 0{y^), 

(24) 
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where 



p ■- 


f^^Po \ 

f^^Po J 


9 


f^qo \ 
f^^qo , 
f^^qo ) 


_ f di2 \ _ f ^22 \ 

■” V 022^22 y ’ y -C22dl2 ) 

with general solution: 


note that we only need wi('d) since the nonlinear 
function in (15) only depends on xi; then by sub¬ 
stituting (34) into (201, we obtain: 


yi= uy2 +gi{yi,y2]p) 
y 2 = -ujyi+ g 2 {yi,y 2 \p) ’ 


(35) 


or 


2/1 = W2/2 + a 2 oyl + aii2/i2/2 + 0022/2 + 0302/1 

+ 0212/12/2 + Oi22/l2/i + 0032/2; 

2/2 = -W 2/1 + ^202/1 + ^112/12/2 + ^022/2 + ^ 302 /i 

+ ^212/i2/2 + ^ 122 / 12/2 + ^032/2- 

(36) 


/i(i?) = e‘=’’’/s: -b M cos(a;i?) -b TV sin(wi?). (30) 


After substituting the general solution into (29) we 
solve for M and N, and then from the boundary 
value problem we solving for K ^ 


C 

ujI 



(31) 


P/i(0) -b Qh{-T) = p - r, 


(32) 


where 


P ■= 


Q:= 



-c, 


(33) 


and r := ( 0 /f 0 /“ 0 /f )^. 

The expressions for wi(0) and wi(—r), necessary 
in (22), are: 

wi(0) = ^ [{Ml + Ki)yl + 2(M3 -b ^ 13 ) 2 / 12/2 


+ (-^5 + Kb)yi 


wi(-r) = i^(e -bMiCOs(a;T) 

- Ni siTL{oJT))yl + 2(e“‘^'^A:|3 

-b M3 cos{ojt) — N3 sm{uJT))yiy2 
+ (e“‘^'^Ar|5 -b M5 cos(clit) 

- A^s sin(a;T))j/| ), 


( 34 ) 


In |3] is computed the coefficient a, which deter¬ 
mines stability of the normal form (36), 


= ^[gr + 9¥+ 91^ + 9?] 


1 

16a; 


9 lH 9 ?+ 9 T)- 9 lH 9 ^^+ 9 r) 


92^ 9f 


■ 2/i°9?° 


(37) 


d'^+3 

where glP = ——— gr(0,0). Periodic orbits near 
oyid3y2 

Hopf bifurcation at the critical eigenvalue A = lo;, 
will be stable if a < 0 and unstable if a > 0. 


4 Numerical Results 

We reproduced some of the computations for the 
Hopf bifurcations in the Fixed point space for 
the case AT > 1 presented in [2], we will com¬ 
pute stability for these bifurcation curves using 
results obtained in the previous section. In fig¬ 
ure (part of figure 10, in [2]) are shown the 
symmetry-preserving bifurcations curves in the pa¬ 
rameter space (/t,T) for K = 1.05, for both cases: 
bifurcations crossing from the left to the right in 
black color, and crossing from the right to the left 
in red color; we also choose three testing point for 
numerical simulation A = (0.3, 6.34), B = (0.3,11), 
and C = (0.4,8.204). 

In figure is shown the coefficient a computed 
using equation ( [37| ), for K = 1.05 in the param¬ 
eter space (/i,r) related to the Hopf bifurcations 
shown in figure the black curve correspond to 
stability of Hopf bifurcations crossing from the left 
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to the right (black curves in figure [^ , as we can 
see, periodic solutions emerging at these Hopf bi¬ 
furcations are all stable (a < 0); the red curve cor¬ 
respond to stability of periodic orbits near Hopf 
bifurcations crossing back from the right to the left 
(red curves in figure [^, they are all unstable for 
« 0.382, and stable for ^ > /i*. 



Figure 1: Symmetry-preserving bifurcations curves 
in Fix(S 7 v) for K = 1.05. In black, bifurcations 
crossing from the left to the right, in red, bifurca¬ 
tions crossing from the right to the left. 


cation crossing from the left to the right in figure, 
which is stable (a < 0), is shown in figure]^ peri¬ 
odic solution near Hopf bifurcations crossing from 
the right to the left at point B, n = 0.3, r = 11, 
which is unstable (a > 0), is shown in figure]^ and 
periodic solution at point C, = 0.4, t = 8.204, 
which is stable (a < 0), is shown in figure all 
initial conditions were set 2 /i( 0 ) = j/ 2 ( 0 ) = 10^. 


X10-* 



0 0.5 1 1.5 2 2.5 3 3.5 4 4.5 


t 

Figure 3: yi{t) at point A, for fi = 0.3 and r = 6.34, 
c.i. 2 /i( 0 )= 2 / 2 ( 0 ) = 10-5. 



Figure 2: Coefficient a computed using eq.p7|, de¬ 
termining stability of Hopf bifurcations in Fix(Sjv) 
for K = 1.05, see figure]^ 


In order to confirm our results, numerical simu¬ 
lations computing yi and ?/2 in equation (36) were 
run for point A, B, and C, using ODE4-5 Mat- 
lab rutine, with time spam 5 x 10^, and maximum 
step size 0.05. Periodic solution yi{t), for point A, 
fj, = 0.3 and r = 6.34, corresponding to Hopf bifur- 





0 0.5 1 1.5 2 2.5 3 3.5 4 


t xio" 

Figure 4: yi{t) at point B, for y, = 0.3 and t = 11, 
c.i. 2 /i( 0 )= 2 / 2 ( 0 ) = 10-5. 


5 Conclusions 

The reduction of the inifinite-dimensional space 
onto the center manifold in normal form, was ap¬ 
plied to the Fixed point space for the Full-phase 
model in order to analyse the stability of simple 
Hopf bifurcations, in both cases, for bifurcations 
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0.5 1 1.5 2 2.5 3 3.5 4 4.5 5 

t xIO^ 


Figure 5: yi{t) at point C, for = 0.4 and t = 
8.204, c.i. j/i(0)=y2(0) = 10-5. 


crossing from the left to the right and for bifurca¬ 
tions crossing in the other way round. We found 
that near all bifurcations that cross from the stable 
region into the unstable region can emerge periodic 
orbits which are stable (a < 0), and, on the other 
hand, we observed that around all bifurcations 
coming back from the right to the left, unstable 
(a > 0) periodic orbits can emerge for y < y,*{K), 
and stable periodic orbits for y > y*{K). 

Although, we computed the coefficient a for a 
specific value of K, the procedure shown in this 
work is valid for all the parameter space where sim¬ 
ple Hopf bifurcations appear. 

Finally, it is important to spotlight some points 
for further research: First, what is the nature of the 
solutions at the special point y = y*{K), at which 
the coefficient a changes sign. Second, analyze sta¬ 
bility of the degenerate Hopf bifurcations at the 
Fixed point space for K = 1, which are codimen¬ 
sion 2, pure imaginary eigenvalue and zero eigen¬ 
value, and third, the stability of the symmetry¬ 
breaking degenerate Hopf bifurcations which have 
multiplicity — 1. 
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